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The competition between collective quantum phases in materials with strongly 
correlated electrons depends sensitively on the dimensionality of the electron 
system, which is difficult to control by standard solid-state chemistry. We have 
fabricated superlattices of the paramagnetic metal LaNi0 3 and the wide-gap 
insulator LaA10 3 with atomically precise layer sequences. Using optical el- 
lipsometry and low-energy muon spin rotation, superlattices with LaNi0 3 as 
thin as two unit cells are shown to undergo a sequence of collective metal- 
insulator and antiferromagnetic transitions as a function of decreasing tem- 
perature, whereas samples with thicker LaNi0 3 layers remain metallic and 
paramagnetic at all temperatures. Metal-oxide superlattices thus allow con- 
trol of the dimensionality and collective phase behavior of correlated-electron 
systems. 
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Since the discovery of high-temperature superconductivity two decades ago, much effort 
has been undertaken to explore and understand the quantum physics of strongly correlated 
electrons in transition metal oxides (TMOs) (i). The electronic phases can exhibit radically 
different physical properties, and a new generation of electronic devices will become possible 
if the competition between these phases can be systematically controlled (2). However, the 
control options offered by conventional solid-state chemistry are limited. The charge carrier 
concentration in a TMO compound, for instance, can be modified by chemical substitution (3), 
but only at the expense of altering the local lattice structure and electronic energy levels in 
an uncontrolled manner. The dimensionality of the electron system, D, is another key control 
parameter, because low-dimensional metals are known to be more susceptible to collective or- 
dering phenomena (including spin- and charge-ordering instabilities as well as unconventional 
superconductivity) than their higher-dimensional counterparts. Some level of dimensionality 
control has been achieved by synthesizing compounds in the Ruddlesden-Popper series of per- 
ovskite structures, which comprise N consecutive TMO layers per unit cell. In principle, the 
dimensionality of the electron system in these materials can thus be tuned from D = 2 to 3 by 
increasing N. In practice, however, the synthesis requirements become rapidly more demanding 
for large N, and many Ruddlesden-Popper phases have turned out to be unstable. 

Recent advances in the synthesis of TMO heterostructures with atomically sharp interfaces 
indicate an alternative route towards control of correlated-electron systems (2). In principle, the 
carrier concentration in a heterostructure can be tuned by a gate voltage in a field-effect arrange- 
ment, without introducing substitutional disorder, and the dimensionality can be modified by 
means of the deposition sequence of electronically active and inactive TMO layers. In practice, 
however, attempts to implement this approach have faced many of the same difficulties encoun- 
tered in the chemical synthesis of bulk materials. For instance, defects created by interdiffusion 
or strain relaxation can influence the transport properties of the interfacial electron system in 
an uncontrolled manner. These difficulties are compounded by the paucity of experimental 
methods capable of probing the collective phase behavior of electrons in TMO heterostruc- 



2 



tures. Whereas ferromagnetism and ferroelectricity can be detected based on the macroscopic 
magnetic or electric field distribution, the identification of two of the most common collective 
ordering phenomena of correlated electrons, namely charge order and antiferromagnetism, in 
TMO heterostructures and superlattices is much more difficult. 

Motivated by the desire to overcome these difficulties and to realize the potential of TMO 
heterostructures in controlling collective quantum phases, we have carried out a comprehensive 
experimental study of superlattices based on the correlated metal LaNi0 3 in which the dimen- 
sionality of the electron system was used as a control parameter, but the influence of epitaxial 
strain and defects was carefully monitored. An extensive body of prior work on bulk nickelates 
provides an excellent background for our study. While bulk LaNi0 3 is a 3D Fermi liquid (4) 
that remains paramagnetic and metallic at all temperatures, other lanthanide nickelates i?Ni0 3 
with smaller electronic bandwidths exhibit collective metal-insulator transitions with decreasing 
temperature (3). In the insulating low-temperature phase, they exhibit a periodic superstructure 
of the valence-electron charge and a non-collinear antiferromagnetic ordering pattern of spins 
on the Ni atoms (5-8). This implies that the itinerant conduction electrons of LaNi0 3 are highly 
correlated and on the verge of localization. Experiments on a controlled number of atomically 
thin LaNi0 3 layers separated by the electronically inactive wide-gap insulator LaA10 3 are thus 
well suited for attempts to control the phase behavior of a correlated-electron system via its 
dimensionality. We have used wide-band spectroscopic ellipsometry to accurately determine 
the dynamical electrical conductivity and permittivity, which (in contrast to the dc conductiv- 
ity) are not influenced by misfit dislocations. Low-energy muons, which are stopped in the SL 
before they reach the substrate, served as a sensitive probe of the internal magnetic field distri- 
bution. Two consecutive, sharp phase transitions in the charge and spin sector revealed by this 
experimental approach demonstrate that the electronic properties of our SLs are determined by 
electron correlations, and not by interfacial disorder. By changing the LaNi0 3 layer thickness, 
we demonstrate full dimensionality control over the collective phase behavior. 

The superlattices (SLs) were grown by pulsed-laser deposition (9, 10) and comprised N 
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consecutive layers of LaNi0 3 and the LaA10 3 . In order to discriminate between the influence 
of dimensionality and epitaxial strain, we have grown SLs on both SrTi0 3 , which induces 
tensile strain in the overlayer, and LaSrA10 4 , which induces compressive strain. Figure 1 shows 
contour maps of the diffracted X-ray intensity distribution in the vicinity of the 103 perovskite 
Bragg peak for three representative samples: N = 4 and N = 2 SLs grown on LaSrA10 4 (001), 
and an N = 2 SL on SrTi0 3 (001). Both the position and the shape of the overlayer reflection 
are strongly affected by inversion of the type of substrate-induced strain (Figs. IB and 1C), 
but remain essentially unchanged by varying the individual layer thicknesses N (Figs. 1A and 
IB). A detailed analysis of the substrate-induced strain and relaxation effects is provided in the 
Online Supplement (10). In the following we show that the transport and magnetic properties 
of the SLs are only weakly influenced by the strain-induced local structural distortions and 
interfacial defects, but qualitatively transformed by varying the number of consecutive unit 
cells within the LaNi0 3 layers. 

The charge transport properties of the SLs were determined by spectral ellipsometry, which 
yields the frequency-dependent complex dielectric function, e(u) = £i(o>) + ie2,{oS), related 
to the optical conductivity a(u) by e(u) = 1 + A7ria(uj)/cj. This method is very sensitive to 
thin-film properties due to the oblique incidence of light, and insensitive to the influence of 
strain-induced extended defects on the current flow through the atomically thin layers (10,11). 
Figures 2A and 2B show the infrared spectra of 62(0^) for N = 4 and 2 SLs grown on LaSrA10 4 
and SrTi0 3 , respectively, which are representative of the in-plane dielectric response of the 
metallic LaNi0 3 layers. The insets show the corresponding temperature dependencies of Ei at 
a fixed photon energy ?vjj = 30 meV. The gradual evolution of e 2 with temperature over the 
far-infrared range confirms that the N = 4 SLs remain metallic at all temperatures. The N = 2 
SLs, on the other hand, show clear evidence of a metal-insulator transition upon cooling, with a 
sharp onset at T MI = 150 K and 100 K for SLs grown on LaSrA10 4 and SrTi0 3 , respectively. 
For T > T M i, the infrared e(u) spectrum of N = 2 SL is well described by a broad Drude 
response e(u) = — w^/(u; 2 + iwy) with a ratio of scattering rate and plasma frequency 
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7/cUpz ~ 0.2 (lower shaded line in Fig. 2A) that is typical for bulk complex oxides. The effective 

mass enhancement m*/m is estimated from the plasma frequency as m* / m = iz^n^ where 

mu) pi 

n-Ni = | x 1-7 x 10 22 cm -3 , by assuming one electron per Ni atom. With cu p i ps 1.1 eV 
[Fig. 2A and Fig. S7 in the SOM (10)] we obtain m* jm 10 which is in good agreement 
with the value for bulk LaNi0 3 obtained from specific heat measurements (12). Using the 
Fermi energy E F = 0.5 eV derived from the thermopower of LaNi0 3 (12) and 7 from the 
Drude model fit to the infrared spectra, we estimate the mean free path as I = ^j2Ep]rrt . 
For the N = 2 and iV = 4 SLs on both substrates we obtain / = 5 - 6 A and 10 - 12 A, 
respectively (10). Remarkably, the mean free path correlates with the individual LaNi0 3 layer 
thickness, testifying to the high quality of the interfaces. 

The charge-carrier localization at lower temperature can be readily identified through a rapid 
drop in e 2 (T) and progressive deviation of €2(00) from the Drude function due to the formation 
of a charge gap. The temperature evolution of the real part of the dielectric function provides 
complementary information about the optical spectral- weight redistribution at T MI . Figures 
2C and 2D show the temperature dependence of the as-measured permittivity, e\, at an energy 
above the gap (?vjj = 0.8 eV). In the metallic phase, e\ decreases with decreasing temperature, 
following the temperature dependence of the scattering rate 7(T). This is characteristic of a 
narrowing of the Drude peak where the spectral weight is removed from the high-energy tail 
and transferred to the far-infrared range near the origin. The charge-gap formation below T M i in 
N = 2 SLs leads to the reverse spectral-weight transfer from the inner-gap region to excitations 
across the gap, and as a consequence, to an increase in e\. 

The consistent temperature evolution of e\ and e 2 over a broad range of photon energies 
demonstrates the intrinsic nature of the charge-localization transition observed in SLs with 
N = 2 and provides clues to its origin . The spectral- weight reduction within the gap can 
be quantified in terms of the effective number of charge carriers per Ni atom and extracted 
from a sum-rule analysis as ASW = w ^ J Q G [<Ji(T m T M i,u>) - (?i(T « T M i,w)]dw, 
where m is the free-electron mass and n Ni the density of Ni atoms. The upper integration limit, 
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Q G ps 0.43 eV, is a measure of the charge gap and can be identified with the equal-absorption 
(or isosbestic) point, where ai(u) curves at different temperatures intersect (Fig. 2E). The 
charge-gap formation on this energy scale can be attributed to a charge-ordering instability, as 
in the case of bulk lanthanide nickelates i?Ni0 3 with smaller i?-ion radius (5-8). The spectral 
weight ASW ps 0.03 below fi G ps 0.43 eV determined for the N = 2 SLs (Fig. 2E) is of the 
same order as, albeit somewhat lower than, the corresponding quantity ASW = 0.058 below 
Qq ~ 0.3 eV reported at the metal-insulator transition in bulk NdNi0 3 , which is known to be 
due to charge order (13). To highlight the analogy to the behavior in bulk nickelates, Fig. 2F 
shows reference measurements on single 100 nm thick films of NdNi0 3 and LaNi0 3 , measured 
under the same conditions as in Figs. 2C and 2D. Since e\ at 0.8 eV for the single NdNi0 3 film 
displays closely similar temperature dependence as found for N = 2 SLs, we conclude that the 
gap formation in the latter case also reflects charge ordering. In NdNi0 3 the metal-insulator 
transition occurs as a first-order transition with a concomitant non-collinear antiferromagnetic 
ordering at T N = T MI (6-8). The thermal hysteresis in the e\(T) curve in Fig. 2F is consistent 
with the first-order character of the transition, with uniform and charge-ordered phases coexist- 
ing over a broad temperature range. In contrast, there is no discernible hysteresis observed in 
e\(T) of iV = 2 SLs (Figs. 2C and 2D), which suggests a second-order transition. Continuing 
the analogy with the bulk nickelate series, one would then expect another second-order transi- 
tion due to the onset of antiferromagnetic ordering at T N < T M i in the iV = 2 SLs, as in i?Ni0 3 
with small R (Lu through Sm). 

In order to test this hypothesis, we carried out low-energy muon spin rotation (LE-/iSR) 
measurements using the fiE4 beamline at the Paul Scherrer Institute (10, 14), where positive 
muons with extremely reduced velocity can be implanted into specimens and brought to rest 
between the substrate and the LaA10 3 capping layer. Since the muons decay into positrons 
preferentially along the spin direction, they act as highly sensitive local magnetic probes. Figure 
3A shows muon decay asymmetry data from a SL with N = 2 at selected temperatures with 
no external field. At T > 50 K, the asymmetry is described by a Gaussian with relatively slow 
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relaxation, a, given by A(0) x exp(— a 2 t 2 /2) (solid lines in Fig. 3A), typical of dipolar magnetic 
fields generated by nuclear moments of La and Al. As the temperature decreases, there is a 
gradual increase in a from 0.17 /is -1 at 250 K to 0.27 fis^ 1 at 20 K. Below 50 K the asymmetry 
drops sharply, and the /iSR spectra can be fitted well by introducing an additional exponential 
relaxation exp(— At). The fast depolarization rate A reaches a value of ps 17//S -1 at 5 K, 
implying a resulting Lorentzian distribution of local fields with half-width-at-half-maximum 
AB = 0.75A/7 M ps 150 G, where 7^ = 2ir x 13.55 MHz/kG is the muon's gyromagnetic 
ratio (10). The fast increase in A with decreasing temperature below 50 K is similar to the 
behavior in bulk NdNi0 3 (15) and (Y, Lu)Ni0 3 (16) below T N , caused by static internal fields 
from ordered Ni magnetic moments. The wide field distribution AB and the absence of a 
unique muon precession frequency reflects the SL structure with several inequivalent muon 
stopping sites in the alternating magnetic (LaNi0 3 ) and nonmagnetic (LaA10 3 ) layers, probably 
compounded by a complex non-collinear spin structure as in the bulk nickelates (15, 16). 

We used 100 G transverse field (TF) measurements to determine the fraction of muons, f m , 
experiencing static local magnetic fields B\ oc > _B TF (i.e. showing no detectable precession 
with uo = 7 M _B TF ) (10). Figure 3C indicates that the N = 2 SL shows a transition from an 
entirely paramagnetic muon environment (f m = 0) to a nearly full volume of static internal 
fields, with a sharp onset at T N ~ 50 K. The magnetic state at 5 K is robust against externally 
applied transverse fields up to 3 kG (not shown), that is the limit of time resolution of our setup. 
The continuous temperature dependence of f m (Fig. 4) and the absence of thermal hysteresis 
indicate that the magnetic transition for the N = 2 SLs is second-order. At the same time, Figs. 
3B and 3D show that SLs with thicker LaNi0 3 layers remain paramagnetic down to the lowest 
temperatures, as in bulk LaNi0 3 . An additional slow exponential relaxation with A = 0.9 fis' 1 
is seen only at T = 5 K (black symbols and curve in Fig. 3B). This results in a small increase 
in relaxation rate, but no loss in asymmetry of the TF /iSR signal (Fig. 3D). The effect is likely 
due to weak dynamical spin correlations that are quenched already in a field of 100 G, in clear 
contrast to the long-range static magnetic order observed in the N = 2 samples. 
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As a local probe, /iSR does not allow definite conclusions about the magnetic ordering 
pattern in the iV = 2 SLs. However, we can rule out ferromagnetism based on an estimate of 
the ordered moment, fi Ni , on the Ni sites from the distribution of local fields experienced by 
the muons. The highest local field at the shortest /i-Ni distance, c/4 (or c*/2 in pseudo-cubic 
notation) 1.92 A (75), is 4-5 times AB, which corresponds to fi Ni > 0.5/i£. If these moments 
were co-aligned in the ordered state, the corresponding total moment M = pLmnrnVsL ^ 
7.7 x 1CT 4 emu would have been readily detected in magnetization measurements. The absence 
of such an effect is confirmed in magnetometric measurements with sensitivity ~ 1CT 7 emu. 

Figure 4 summarizes the phase behavior of the SLs with N = 2, which undergo a sequence 
of two sharp, collective electronic phase transitions upon cooling. We have provided strong 
evidence that the two transitions correspond to the onset of charge and spin order. By showing 
that the N = A counterparts remain uniformly metallic and paramagnetic at all temperatures, 
we have demonstrated full dimensionality control of these collective instabilities. The higher 
propensity towards charge and spin order in the two-dimensional systems probably reflects 
enhanced nesting of the LaNi0 3 Fermi surface. The phase behavior is qualitatively similar 
to the one observed in bulk i?Ni0 3 with small radius of the R anions, which results from 
bandwidth narrowing due to rotation of Ni0 6 octahedra, but the transition temperatures and 
the order parameters are substantially lower, probably because of the reduced dimensionality. 
Since the transitions occur in the N = 2 SLs irrespective of whether the substrate-induced 
strain is compressive (Fig. IB) or tensile (Fig. 1C), structural parameters such as rotation and 
elongation of the NiC>6 octahedra can be ruled out as primary driving forces. We note, however, 
that the infrared conductivity is higher (Figs. 2A and 2B) and the transition temperatures are 
lower (Fig. 4) in the N = 2 SL grown under tensile strain. The more metallic response of 
these SLs, compared to those grown under compressive strain, may reflect a widening of the Ni 
3d bandwidth and/or an enhanced occupation of the Ni d x 2_ y 2 orbital polarized parallel to the 
LaNiC>3 layers. A small orbital polarization was indeed detected by soft x-ray reflectometry in 
our superlattices (17). This indicates further opportunities for orbital control of the collective 
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phase behavior of the nickelates, which may enable experimental tests of theories predicting 
high-temperature superconductivity (18, 19) or multiferroicity (20) in these systems. 
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Figure 1: Reciprocal-space maps of 100-nm-thick LaNi0 3 (N u.c.)|LaA10 3 (N u.c.) super- 
lattices grown under compressive strain on LaSrA104 (001) with (A) N = A and (B) N = 2 
and (C) under tensile strain on SrTi0 3 (001) with N = 2. The black vertical lines indicate the 
in-plane (Q x ) position of the LaSrA10 4 (109) and SrTi0 3 (103) reflections. The strain state 
of the perovskite epilayers is identified by the intensity distribution in the vicinity of the (103) 
layer Bragg peak and its superlattice satellite, which are delineated by solid- and dashed-line 
triangles, respectively. The reciprocal spacings of 103 strain-free pseudo-cubic LaNi0 3 and 
LaA10 3 are indicated by the red circles. The red arrows point towards the origin. 
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Figure 2: (A,B) 62(00) spectra of the N = 2 (solid lines, circles) and N = 4 (dashed lines, 
triangles) SLs on (A) LaSrA10 4 and (B) SrTi0 3 substrates measured at representative temper- 
atures. The shaded lines in (A) represent the Drude model fit to e 2 (co) at 175 K for the N = 2 
(N = 4) SL with oo p i = 1.05 eV (1.10 eV) and 7 = 200 meV (90 meV). The insets provide 
the corresponding temperature dependencies of £2 at a photon energy of hw = 30 meV for the 
N = 2 (circles) and N = 4 (triangles) SLs. (C,D) Temperature dependence of the as-measured 
pseudo-dielectric permittivity e\ at hu = 0.8 eV in the N = 2 (blue) and N = 4 (black) SLs 
on (C) LaSrA10 4 and (D) SrTi0 3 . (E) The difference between the optical conductivity spectra 
<Ti(100 K, u) and o"i(10 K, of) (shaded area) quantifies the reduction of the effective charge 
density, ASW ~ 0.03 per Ni atom, within the gap energy range below 0.43 eV at the charge 
ordering transition in the N = 2 superlattice on SrTi0 3 . (F) Temperature dependence of e\ at 
0.8 eV of reference 100 nm films of LaNi0 3 (black) and NdNi0 3 (blue). 
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Figure 3: The time evolution of the zero-field muon spin polarization at various temperatures 
for the (A) N = 2 and (B) N = 4 superlattice on LaSrA10 4 . (C) and (D) Muon spin relaxation 
spectra in a weak transverse magnetic field of 100 G in the superlattices as in (A) and (B), 
respectively. (A) to (D) use the color coding in legend (B). 
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Figure 4: Temperature dependencies of the fraction of muons experiencing static local magnetic 
fields, f m , and the normalized permittivity e* at 0.8 eV in the N = 2 superlattices. The black 
and blue arrows mark the magnetic (T N ) and metal-insulator (T MI ) transition temperatures for 
the superlattices on LaSrA10 4 and SrTi0 3 , respectively. 
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Methods and materials 

Sample preparation and characterization 

High-quality superlattices (SLs) composed of N u.c. thick consecutive layers of LaNi0 3 and 
LaA10 3 were grown by pulsed-laser deposition from stoichiometric targets using a KrF excimer 
laser with 2 Hz pulse rate and 1.6 J/cm 2 energy density. Both compounds were deposited in 0.5 
mbar oxygen atmosphere at 730°C and subsequently annealed in 1 bar oxygen atmosphere at 
690°C for 30 min. We have grown SLs on two kinds of single-crystalline substrates: SrTi0 3 , 
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Table SI: Average lattice constants of 100 nm thick N = 2 SLs grown on (001) -oriented SrTi0 3 
and LaSrA10 4 substrates (determined from the main (103) layer Bragg peak positions in Figs. 
IB and 1C) in comparison with the lattice constants of strain-free pseudo-cubic LaNi0 3 and 
LaA10 3 and the same substrates. 





SL on LaSrA10 4 


SL on SrTi0 3 


LaNi0 3 


LaA10 3 


LaSrA10 4 


SrTi0 3 


a, b (A) 


3.769 


3.845 


3.837 


3.789 


3.756 


3.905 


c(A) 


3.853 


3.790 


3.837 


3.789 


12.636 


3.905 



which induces tensile strain in the overlayer, and LaSrA104, which induces compressive strain 
(see Table SI). All substrates were lOmmx 10mmx0.5 mm or 5 mmx5 mmx0.5 mm (001)- 
oriented plates with a miscut angle < 0.1°. We chose to work on 100 nm thick SLs in order 
to enhance the dielectric response and to confine the muon stopping distribution within the SL. 
The chosen thickness range also allows us to avoid complications arising from initial growth of 
TMO layers on a substrate (SI). The growth rates for the individual layers were controlled by 
counting laser pulses in combination with feedback from high-resolution x-ray diffraction mea- 
surements. The crystallinity, superlattice structure, and sharpness of the interfaces (with rough- 
ness < 1 u.c.) were verified by momentum-dependent x-ray reflectivity and high-resolution hard 
x-ray diffraction scans which revealed, besides the perovskite Bragg reflections, satellite peaks 
due to the long-range multilayer superstructure and Kiessig fringes caused by total-thickness 
interference. 

Representative scans along the specular truncation rod are shown in Fig. SI for samples 
grown on the different substrates with different individual layer thicknesses N (u.c), and total 
thicknesses D (A). Symmetrically around the (001) layer Bragg peak one can see superlattice 
satellites and M — 2 thickness fringe maxima, where M is the number of bilayer repetitions. 
The position of the satellites corresponds to a LaNi0 3 (N u.c.) |LaA10 3 (N u.c.) bilayer thick- 
ness of 30 ± 1 A and 15.5 ± 0.5 A for the N = 4 (Figs. S1A and SIB) and N = 2 (Figs. 1C 
and ID), respectively, that is in a good agreement with the 2A^c value, where c is the average 
epilayer lattice constant in Table SI. Accordingly, the Kiessig fringes in Figs. SI and S2 cor- 
respond to the total thickness M x 2Nc. The thickness fringes for the 100 nm thick N = 2 
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Figure SI: High-resolution x-ray diffraction measured with lOkeV synchrotron radiation at 
the MPI-MF beamline of the ANKA facility at the Karlsruhe Institute of Technology for the 
(A, B) 30 nm (N = 4) and (C, D) 100 nm (N = 2) thick superlattices on (A,C) SrTi0 3 and 
(B, D) LaSrA10 4 . The thickness of the SL in (D) is determined from the hard x-ray reflectivity 
measurements in Fig. 2D. The 100 nm thick samples were used for low-energy muon spin 
rotation and ellipsometry experiments. 



SL on LaSA10 4 are damped at higher I values in Fig. SID, but well resolved in Fig. S2D. 
The x-ray reflectivity shown in Fig. S2 was also used to characterize the superlattice structure 
and sharpness of the interfaces. From fits to the reflectivity, using the Parratt algorithm and 
tabulated values for the optical constants (S2,S3, S4), we obtained the thickness {d^Ao, diNo) 
and roughness {o^ao-, <Jlno) of the individual layers. Using a minimal set of fitting parameters 
(assuming M identical LaNi0 3 and LaA10 3 layers), we show in Fig. S2 a good description 
of the data. The roughness parameters are all around 1 u.c or less and represent values aver- 
aged over a large area of ~(10x 1) mm which corresponds to the x-ray spot size and inevitably 
contains planar defects such as stacking faults. This indicates the presence of atomically flat 
and abrupt interfaces. Some of the samples were also checked by high-resolution transmission 
electron microscopy (TEM), providing a local picture of the atomic stacking sequence. In Fig. 
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Figure S2: Hard x-ray reflectivity measured with Cu Ka radiation and fits for the (A) 23 nm 
and (B) 100 nm thick N = 2 superlattices on SrTi0 3 and 100 nm thick (C) iV = 2 and (D) 
N = A superlattices on LaSrA104. The samples in (C-D) were used for low-energy muon spin 
rotation and ellipsometry experiments. 

S3 a high-angle annular dark-field image of a LaNiC>3 (2 u.c.) |LaA103 (1 u.c.) SL is shown. In 
this imaging mode, also known as Z-contrast, the contrast is proportional to Z n , where Z is the 
atomic number and n is about 1.7. Subsequent dark (marked by arrows) and bright layers show 
a chemical variation of the layer system. In this example a sequence of two LaNiOs layers and 
one LaA10 3 layer is visible which shows that even single layers can be deposited without dis- 
tinct intermixing. The superior quality of our samples is also supported by resonant reflectivity 
measurements performed on a sample grown under the same conditions. The analysis of those 
data allowed some of us to determine the atomic-layer resolved orbital polarization in these 
superlattices (S5). 



Substrate-induced strain and relaxation effects 

In general, the physical properties of thin films are strongly influenced by substrate-induced 



Figure S3: High-angle annu- 

lar dark field image of the 
LaNi0 3 (2 u.c.)|LaA10 3 (1 u.c.) 
superlattice. Subsequent dark (marked 
by arrows) and bright layers show the 
chemical variation of the layer system. 

strain and relaxation effects. It has thus far proven difficult to separate the influence of the di- 
mensionality from that of other parameters such as the strain-induced local structural distortions 
and interfacial defects. In order to discriminate between these effects we chose to work on SLs 
grown on both SrTi0 3 , which induces tensile strain in the overlayer, and LaSrA10 4 , which in- 
duces compressive strain. Our comprehensive reciprocal- space mapping (RSM) measurements 
(S6) supplemented by high-resolution TEM micrographs verified that strain and relaxation ef- 
fects are strongly affected by inversion of the type of substrate-induced strain, but remain essen- 
tially unchanged by varying the individual layer thicknesses. In our study, we show that, on the 
contrary, the transport and magnetic properties of the SLs are almost unaffected by inversion 
of the type of substrate-induced strain, but qualitatively transformed by varying the number of 
consecutive unit cells within the LaNi0 3 layers. Since the metal-insulator and spin-ordering 
transitions occur in the N = 2 SLs irrespective of whether the substrate-induced strain is com- 
pressive or tensile, strain-induced local structural distortions and interfacial defects are ruled 
out as primary driving forces. 

Figure 1 of the main text shows contour maps of the diffracted X-ray intensity distribu- 
tion in the vicinity of the 103 perovskite Bragg peak for three representative samples: N = 4 
and iV = 2 SLs grown on LaSrA10 4 (001), and an N = 2 SL on SrTi0 3 (001). The anal- 
ysis of the averaged in-plane and out-of-plane lattice constants (Table SI) indicates that com- 



Figure S4: Reciprocal- space maps in the vicinity of the symmetric (004) peak of the 100 nm 
thick N = 2 superlattices on (A) SrTiC>3 and (B) LaSrAlC>4 substrates. The relaxation triangle 
is highlighted with a red line in (A). The angle (3^2° quantifies the amount of gradual relax- 
ation the SL has. (C) Horizontal cuts along the indicated in (B) q z values. The separation of the 
two twin peaks reveal the formation of twin domains. 

pressive strain reduces the in-plane lattice parameter by Aa/a « 1.8 % relative to the bulk 
LaNiC>3 lattice, whereas tensile strain results in a reduction of the out-of-plane lattice constant 
by Ac/c Rj 1.2 %. These two types of local distortions in the perovskite structure are accommo- 
dated by rotations of the Ni0 6 octahedra about different Cartesian axes (S7), which, in turn, ex- 
ert an inequivalent influence on the LaNiC>3 electronic structure. A distribution of the diffracted 
intensity near the epilayer reflection for SLs grown on LaSrA10 4 has a characteristic triangular 
shape, with dispersion along the in-plane (Q x ) direction towards the 103 Bragg reflection of 
strain-free bulk LaNiC>3. This is in contrast to the tensile-strained SLs grown on SrTiC>3, where 
the strain relaxation is characterized by nearly elliptical contour lines close to the 103 Bragg re- 
flection of cubic LaAlOs. The tensile strain of SrTiC>3 is (67 ± 3)% relaxed, and by comparing 
with SLs of less total thickness, we identified a faint gradient-profile-relaxation effect as func- 
tion of overlayer thickness, similar to the behavior observed in semiconductor heterostructures 
(S8, S9). In addition to Fig. 1 of the main text, Fig. S4 shows RSMs of the symmetric 004 per- 
ovskite Bragg peak measured with synchrotron radiation at the MPI-MF beamline of the ANKA 
facility at the Karlsruhe Institute of Technology. The diffracted x-ray intensity distribution for 
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Figure S5: High-resolution TEM micrographs of LaA10s|LaNi03 SLs on (A) LaSrA104 and 
(B) SrTiOs substrates. Defects are marked by arrows. The inset in (A) shows a magnified area 
close to a planar defect. 

the N = 2 SL on SrTiOs (Fig. S4A) exhibits the triangular shape described in Refs. S8 and 
S9. The effect of triangular relaxation was not observed in the thinner SL (not shown). Because 
the lattice constants of the thin SL are almost equal to the ones of the thick sample, the SLs 
grown on SrTiOs seem to relax abruptly at the beginning of the growth. Further, from the peak 
shape evolution, the SL gradually relaxes the tensile strain. The subtle thickness evolution of 
the layer's relaxation indicates that it is the substrate surface where abrupt strain- adapting mech- 
anisms take place. The effect of tensile strain on TMO heterostructures may produce oxygen 
vacancies (S10), which give rise to a different valence state of the Ni ion at the substrate inter- 
face (SI). Figures S4B and S4C show that the distribution of the diffracted intensity near the 
epilayer reflection for SLs grown on LaSrA10 4 has a double-peak splitting along the in-plane 
(Q x ) direction. This intensity pattern (only seen in thicker SLs grown on LaSrA10 4 ) suggests 
the formation of twinning domains, as described in Refs. Sll, S12. The two different relax- 
ation mechanisms in the perovskite structure are confirmed by TEM measurements performed 
on samples grown under the same conditions as in our study. Figure S5 shows high-resolution 
TEM micrographs (recorded by a JEOL JEM4000FX microscope) of the LaNi0 3 - LaA10 3 
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layer systems. In the case of the LaSrA10 4 substrate (Fig. S5A) planar defects are visible 
(marked by arrows) which are oriented perpendicular to the substrate plane and extend through 
the entire SL. As shown in the magnified inset image, the stacking sequence changes at these 
faults (yellow broken line). The size of the defect- free blocks varies between 15 and 50 nm. 
The micro structure of the layer system on the SrTi0 3 substrate (Fig. S5B) only very occasion- 
ally shows planar defects. Instead, localized defects are found close to the substrate (marked by 
arrows). These defects can be associated with the creation of oxygen vacancies and changes in 
the oxygen coordination of Ni ions at the substrate interface. Recent photon energy-dependent 
hard x-ray photoelectron spectroscopy measurements on some of our samples have confirmed 
that the initial growth on the SrTi0 3 surface leads to the Ni 2+ valence state (S13). The oxygen 
vacancy formation energy gradually decreases with increasing the in-plane perovskite lattice 
spacing (S14), which can explain the marked difference in the oxygen vacancy concentration 
in thin films grown under tensile or compressive strain (S10). Nevertheless, in our study, the 
temperature-induced phase transitions occur in the N = 2 (but not in N = 4) SLs irrespective of 
whether the substrate-induced strain is compressive or tensile, which clearly distinguishes these 
transitions from those in highly oxygen deficient LaNi0 3 _(5 (5 > 0.25) (S15, S16). Moreover, 
the reduced insulating phases require more than 1/3 of divalent Ni 2+ in square planar (vs. per- 
ovskite octahedral) sites. Based on the detailed characterization of our samples by means of 
XRD, XAS, RSM, HAXPES, and TEM we can definitively rule out such a scenario. 

In conclusion, our analysis confirms the excellent quality of the synthesized SLs, which 
exhibit abrupt interfaces and excellent crystallinity. Defect- free, atomically precise 15-50 nm 
blocks are separated by ~ 1 u.c. stacking faults. These planar defects are inevitably caused by 
strain relaxation effects, and can block the current flow through the atomically thin layers. We 
have therefore used advanced local probes, such as spectroscopic ellipsometry and low-energy 
muons, to study the intrinsic electronic transport and magnetic properties of the heterostruc- 
tures. 
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Figure S6: Experimental (open circles) and best-fit calculated (solid lines) ellipsometry spectra 
of the N = 2 and N = 4 SLs on (A,B) LaSrA10 4 at T = 175 K and (C,D) SrTi0 3 at T = 100 
K. The angle of incidence of the polarized light was $j = 82.5°. Ellypsometry spectra of the 
bare substrates measured at $j = 77.5° are shown for comparison (black solid lines). The 
gray shaded area in (C,D) indicates the region where the data analysis is affected by dielectric 
microwave dispersion of the ferroelectric soft mode of SrTiC>3. 



Spectroscopic ellipsometry measurements and data analysis 

We have used wide-band spectroscopic ellipsometry to accurately determine the dynamical 
electrical conductivity and permittivity of the SLs. The distinct advantages of ellipsometry 
are as follows, (i) In contrast to dc transport experiments, this method exposes the intrinsic 
electrodynamic response of the SLs, which is not influenced by the substrate, contacts, and ex- 
tended defects, (ii) As a low-energy spectroscopic tool, it serves to determine critical parame- 
ters of the metal-insulator transition such as the energy gap and the density of carriers localized 
below Tmi- (iii) In comparison with other spectroscopic techniques, ellipsometry yields the 
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Figure S7: Best-fit model functions ei(u) and 82(00) for the N = 2 and N = 4 SLs on (A) 
LaSrA10 4 at T = 175 K and (B) SrTi0 3 at T = 100 K, as obtained by inversion of the 
ellipsometric parameters in Fig. S6. The shaded lines represent the Drude model simultaneous 
fit to both ei(u) and e 2 (co) with parameters cu p i and 7 described in the legends. The gray shaded 
area in (B) indicates the region where the model fitting curves deviate significantly from the 
data. 



frequency-dependent complex dielectric function without the need for reference measurements 
and Kramers-Kronig transformations, (iv) Variable angle ellipsometry is very sensitive to thin- 
film properties due to the oblique incidence of light, and it is generally used to derive optical 
constants of thin films and complex heterostructures (S17). 

The experimental setup comprises three ellipsometers to cover the spectral range of 12 meV 
to 6.5 eV. For the range 12 meV to 1 eV, we used a home-built ellipsometer attached to a 
standard Fast-Fourier- Transform Bruker 66v/S FTIR interferometer. The far-infrared measure- 
ments were performed at the infrared beamline IR1 of the Angstrm Quelle Karlsruhe ANKA 
synchrotron light source at the Karlsruhe Institute of Technology. For the mid-infrared mea- 
surements, we used the conventional glow-bar light source from a Bruker 66v/S FTIR. Finally, 
temperature dependencies of the pseudo-dielectric permittivity el at huj = 0.8 eV were mea- 
sured with a Woollam variable angle spectroscopic ellipsometer (VASE) equipped with an ultra 
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Figure S8: Temperature dependence of the dc (A) resistivity and (B) conductivity of the N = 2 

(green) and N = 4 (blue) SLs on LaSrA10 4 . Solid squares represent the conductivity obtained 

from the Drude model parameters in the legend of Fig. S7A, which exceeds the corresponding 
a dc (T = 175 K) yalues by less than 20 %. 

high-vacuum cold-finger cryostat operated at < 5 x 10~ 9 mbar chamber pressure. 

The inherent capacity of Woollam VASE ellipsometers to measure relative changes of the 
dielectric function on the order of 10~ 2 was boosted to an unprecedented level of 1CT 4 us- 
ing temperature-modulation measurements of the dielectric constant at particular photon en- 
ergies. The ellipsometric angles ^ and A are defined through the complex Fresnel reflection 
coefficients for light polarized parallel (r p ) and perpendicular (r s ) to the plane of incidence, 
tan\l/ e lA = r p /r s . Figure S6 shows representative infrared spectra of ty(u) and A(cu) for 
the = 4 and 2 SLs and for the bare LaSrAlC>4 and SrTiC>3 substrates. The details of the 
data analysis have been discussed elsewhere (S18). The SLs were treated as single-layer films 
according to an effective-medium approximation with a mixture of the nickelate and aluminate 
layers. A wavelength-by-wavelength regression procedure has been employed to extract the 
real and imaginary parts of the dielectric function (SI 9). Figure S7 shows the best-fit model 
functions e\(uS) and ezioS) obtained by inversion of the ellipsometric parameters in Fig. S6. The 
infrared spectra are well described by a broad Drude response e(uj) = — u&J (cu 2 + iuj) 
with a ratio of scattering rate and plasma frequency 7/0^ ~ 0.1 — 0.2 that is typical for bulk 
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complex oxides. The parameters in the Drude fit are well constrained, because both Si(uj) and 
e 2 (uj) are available. The deviation of the Drude fit from the measured Ei(u) and 82(0*1) below 30 
meV in Fig. S7B (gray shadow area) can reflect the uncertainty in the inversion procedure for 
SLs on SrTi0 3 due to the microwave dispersion of the ferroelectric soft mode of SrTi0 3 (S20) 
and/or due to the presence of a dead layer with reduced conductivity at the substrate interface 
(Fig. S5B). This low-energy uncertainty does not, however, affect the relative spectral weight 
reduction, ASW ~ 0.03 (±10%) per Ni atom within the gap energy range below Vic ~ 0.43 
eV, at the metal insulator transition in the iV = 2 SL on SrTi0 3 . 

The effective mass enhancement m*/m is estimated from the plasma frequency as 

*/ 47re 2 n 11.7 

where n — \ x 1.7 x 10 22 cm -3 , by assuming one electron per Ni atom. We note that u p i is 
almost independent of N, implying the volume fraction of the metallic LaNi0 3 layers remains 
the same in all SLs. With oo p i w 1.1 eV, as derived from the Drude model fit in Fig. 2A, 
we obtain m*/m pa 10 which is in good agreement with the value for bulk LaNi0 3 from the 
specific heat measurements (S21). Using the Fermi energy Ep = 0.5 eV derived from the 
fhermopower of LaNi0 3 (S21), we estimate the Fermi velocity as 



2Ep m - 
VF = C \ — - — ps 1.33 x 10 7 cm/s. (2) 
V mc 2 m* 

The mean free path, I, can be estimated from 

l [A] = v F r = J*- = 6.57 x 10-»^J^ « - ™ (3) 
27TC7 7 [meVJ 7 [meVJ 

With 7 p^200 meV (90 meV), as derived from the Drude model fit in Fig. 2A and Fig. S7A, 
we obtain I = 4.4 A (9.7 A) for the iV = 2 (iV = 4) SL on LaSrA10 4 . For the iV = 2 
(N = 4) SL on SrTi0 3 we estimate / = 6.4 A (12 A), respectively. Remarkably, the mean free 
path correlates with the individual LaNi0 3 layer thickness, testifying, along with the constant 
volume fraction of the metallic layers, to the atomic quality of the interfaces . 
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Our results indicate that, even in the iV = 2 samples at T > T MI , the conductivity of 
the LaNi0 3 layers exhibits a clearly metallic temperature and frequency dependence. We de- 
fine T MI as the temperature at which the temperature derivatives of both e 2 (T) (Figs. 2A 
and 2B) and £i(T) (Figs. 2C and 2D) change sign. The consistent temperature evolution 
of e\ and e 2 over a broad range of photon energies demonstrates the intrinsic nature of the 
charge-localization transition observed in SLs with N = 2. In the to — > limit this criterion 
is analogous to a sign change of the temperature derivative of the dc resistivity, dp/dT, ob- 
served at Tmi in bulk _RNi0 3 . This is in contrast to results of recent dc electrical resistivity 
measurements where the insulating behavior of 2 u.c. thick LaNi0 3 is attributed to variable 
range hopping transport (S22) or film-substrate interface effects (S23, S24). For a thickness 
of iV = 2 u.c. the behavior of (LaNi0 3 ) N /(SrMn0 3 )2 SLs is insulating over the entire tem- 
perature range, whereas the N = 4 u.c. SL is metallic with an upturn in resistivity below 
50 K. Even in the iV = 4 u.c. metallic sample, the mean free path I is estimated to be less 
than a single unit cell (S22). This suggests that Anderson localization induced by disorder is 
responsible for the insulating behavior in these systems, in contrast to the sharp temperature 
dependence observed in our SLs that indicates a metal-insulator transition driven by collective 
interactions. Ultrathin single films of LaNi0 3 show a crossover from metallic to insulating 
behavior at a larger thickness (S23, S24), which varies from 6 u.c. to 13 u.c. depending on 
the substrate. We argue that the presence of planar stacking fault defects and a dead layer with 
reduced conductivity at the substrate interface, as discussed above (Fig. S5), makes the analysis 
of the temperature-dependent resistivity measurements challenging and inconclusive about the 
conduction mechanism of ultrathin LaNi0 3 films. Nevertheless, in order to directly compare 
our results with those reported in Refs. S22, S23, and S24, we have also performed dc resis- 
tivity measurements on the iV = 2 and N = A SLs on LaSrA10 4 . Figure S8 shows that the 
dc resistivity of the N = 2 SL exhibits a crossover from metallic to insulating phase behavior 
below T M i ~ 150 K . The sharp temperature dependence in the insulating state does not fit 
to the stretched exponential function (S22) and can not be attributed to variable range hopping 
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transport. Figure S8 also shows that our dc and optical conductivity data (Figs. 2A and S7A) 
are in close agreement. This resemblance once more indicates a low density of stacking faults 
in our samples. 

Low-energy ^SR instrumentation and data analysis 

Low energy muon spin rotation/relaxation (LE-/iSR) uses ~ 100% spin polarized positive 
muons of tunable keV-scale energy to study local magnetic properties of thin films or het- 
erostructures as a function of the muon implantation depth. The details of the data acquisition 
and analysis have been described elsewhere (S25, S26) . More details of the LE-yuSR meth- 
ods and apparatus can be found on the website of the LEM group at Paul Scherrer Institute 
(S27). This technique has been recently successfully applied to the case of magnetic ultra- 
thin films (S28) and wires (S29). Figure S9 shows the muon stopping profile calculated for 
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Figure S9: Muon stopping profile in the o 

N = 2 SL showing the calculated prob- ;g 

ability that /i + with an implantation en- | 
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and 15 keV (red) comes to rest at a cer- 
tain depth near the surface. 

LaA10 3 |LaNi0 3 SLs using the Monte Carlo algorithm TRIM.SP (S30, S31) . 

In our study we found that varying the stopping distribution of ji + on the scale of about 
50-800 A through the control of the implantation energy between 5-15 keV had no effect on 
the fiSR spectra. The experimental LE-^SR curves in Figs. 3A to 3D were measured with 
muons of energy 10 keV, which are implanted at a mean depth of 45 nm. The initial asymmetry, 
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Figure S10: Zero-field fiSR function observed in the iV = 2 SL on LaSrA10 4 at 5 K. The solid 
lines represent the best-fit curves for (A) the two-component model function described in the 
manuscript and (B) Eq.(4), respectively. 

A(0) « 0.18, is smaller than the asymmetry of the LE-/iSR setup of ~ 0.27, because only 2/3 
of the muon beam with a diameter of about 2 cm hit the sample with an area of 1 x 2 cm 2 . The 
sample was surrounded by a Ni-coated sample holder, which causes a very fast depolarization 
(< 0.06 fis) of muons missing the sample. 

The obtained spectra /iSR spectra yield the probability distribution of the local magnetic 
field at the muon sites. As a local probe, fiSR does not allow definite conclusions about the 
magnetic ordering pattern in the N = 2 SLs. However, we rule out ferromagnetism based 
on an estimate of the ordered moment on the Ni sites from the fiSR lineshape, fim > 0.5/zb 
(see the main text of the manuscript). If these moment were co-aligned in the ordered state, 
the corresponding total moment M = fimn-rnVsL > 7.7 x 10~ 4 emu would have been readily 
detected in magnetization measurements. The absence of such an effect, which we confirmed 
in magnetometric measurements with sensitivity ~ 10~ 7 emu. 

We can also rule out a spin-glass state as the ground state of N = 2 SLs, bearing in mind 
that oxygen deficient LaNi0 2 . 75 exhibits spin-glass like behavior at low temperatures (S15). A 
spin-glass state develops gradually due to randomly fluctuating local moments. In this case, 
the spin relaxation function should be exponential with a unique rate already at temperatures 
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above about four times the actual glass transition temperature, ~ 80 - 100 K (S32), which is at 
variance with the sharp temperature onset of the local moment observed in our data (Fig. 3A 
and solid squares in Fig. 4 of the manuscript). A similarly sharp transition was very recently 
observed by x-ray magnetic circular dichroism (XMCD) measurements in a magnetic field of 
5T on a sample with N = 2 grown under the same conditions (S33). 

Additional evidence against a spin glass state can be derived from an analysis of the muon 
relaxation function. At low temperature, the spin-glass relaxation function in zero field can be 
described by (S34) 

A{t) = A [^exp(- V / V) + \{l - - 7 =^=)exp(- v / A^T^] (4) 

with a = y/q<J s and \ d = 4ct^(1 — q)/v, where q is the Edwards- Anderson order parameter 
with the purely static and dynamic limits, q = 1 and q = 0, respectively, a s is the static width 
of local fields at the muon site, and v is the rate of the randomly fluctuating moments. This 
form of the relaxation function is expected for fx + in coexisting static and dynamic random 
local fields. The fit of Eq. (4) to the time evolution of the zero-field muon spin polarization for 
the N = 2 SL on LaSrA10 4 at 5 K (Fig. S10) gives reasonable parameters, i.e. close to the 
static limit with q 0.993, a s = 11 — 15 /is -1 , and v « 2 /is -1 . Nevertheless, the simpler 
two-component model function, as described in the manuscript, provides a better fit to the data 
below 0.2 /is than the spin-glass function of Eq. (2). The analysis is consistent with long-range 
static antiferromagnetic order and confirms the conclusion of our manuscript. 
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